PERFORMANCE  ANALYSIS 
OF  A 

HARTMANN  WAVEFRONT  SENSOR 
USED  FOR  SENSING 

ATMOSPHERIC  TURBULENCE  STATISTICS 


Acknowledgements 


Table  of  Contents 


List  of  Figures 


5 

II 

PO 

2-a 

Lo/D 

5 

II 

Po 

2-a 

a  / 

S  ^ 

II 

Po 

2-a 

a 

Lo 

D 

OO 

00 

s  ^ 

II 

Po 

2-a 

a 

Lr 

D 

00 

s  ^ 

II 

PO 

2-a 

a 

Lo 

D 

00 

s  ^ 

II 

PO 

2-a 

a 

Lo 

D 

00 

s  ^ 

II 

Po 

2-a 

Lo/D 

s  ^ 

Lo 

D 

OO 

00 

Ds  u 

II 

PO 

2-a 

a 

Lo 

D 

oo 

00 

Ds  u 

II 

PO 

2-a 

a 

Lo 

D 

00 

Ds  u 

II 

Po 

2-a 

a 

Lo 

D 

00 

Ds  u 

II 

Po 

2-a 

a 

Lo 

D 

00 

Ds  u 

D/po 

00 


a 

Lo/D 


II 

PO 


a 


a 


Ln 

D 


oo 


n  2— a 

Ds  U  - 
*  Po 

Lq/d 

a 

eo 

^sm  Pi 

-  •  Ai) 

(Jn,VTt,pA 

Ai)  Ai) 

5  nm  5  Pn  5 

Aj) 

Q}^ 

N 

P 

.  z?,  . 

0 

Lo/D 

a 

VT 

•  D, 

b 

0 

N 

p 

.  z?,  . 

0 

Lo/D 

a 

VT 

•  D, 

b 

0 

N 

p 

.  D,  . 

0 

Lo/D 

a 

VT 

• 

b 

0 

N 

p 

.  D,  . 

0 

LolD 

a 

VT 

•  D, 

b 

0 

N 

p 

.  D,  . 

0 

LolD 

a 

VT 

• 

b 

0 

N 

p 

.  D,  . 

0 

Lo/D 

a 

VT 

• 

b 

0 

N 

p 

D,  .  ° 

s 

a 

VT 

•  D, 

.  °  LqJD  00 

N 

p 

.  D,  . 

0 

VT  JD 

a 

V 

°  L^/D 

00 

N 

P 

.  D,  . 

0 

V 

a 

■u 

.  D  Lq/D  00 

r^o) 

S  {a2«2  u  }  /£  {a2«2 

} 

£  {03(13  u  }  /£  {0303 

} 

£  {0303  u  }  /£  {0303 

} 

a 


S 

E 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 


^j^j'(u) 

^j^j'(u) 

^j^j'(u) 

^j^j'(u) 

^j^j'(u) 

^j^j'(u) 

^j^j'(u) 

^j^j'(u) 

^j^j'{u) 

^j^j'{u) 

^j^j'{u) 

^j^j'(u) 

^j^j'{u) 


(2-«) 

PO 

(2-«) 

Po 

(2-«) 

PO 

(2-«) 

Po 

(2-«) 

PO 

(2-«) 

Po 

(2-«) 

PO 

(2-«) 

Po 

(2-«) 

PO 

(2-«) 

PO 

(2-«) 

Po 

(2-«) 

PO 

(2-«) 

Po 

(2-«) 

PO 

(2-«) 

Po 

(2-«) 

PO 

(2-«) 

Po 

(2-«) 

PO 

(2-«) 

Po 

Oo  °  Lo/D 

00  °  Lo/D  .  a 

00  °  Lq/D  oo  a 

00  °  Lo/D  .  a 

00  °  Lo/D  .  a 

00  °  Lo/D  oo  a 

^0  °  Lo/D  .  a 


00  °  Lo/D  .  a 

0Q  °  Lq/D  oo  a 
0Q  °  Lq/D  .  a 

00  °  Lo/D  .  a 

00  °  Lo/D  oo  a 

00  °  Lo/D  .  a 

00  °  Lo/D  .  a 

00  °  Lo/D  oo  a 

00  °  Lo/D  .  a 

00  °  Lo/D  .  a 


00  °  Lo/D  oo  a 

00  °  Lo/D  .  a 


S 

£ 

£  djOijif^^^ 
£  djOijif^^^ 
£  djOijif^^^ 

£  CljClji^j;^^ 
£  CljClji^j;^^ 

£  djOijif^^^ 
£  djOijif^^^ 

£  CljClji^j;^^ 

£  djOijif^^^ 
£  djOijif^^^ 
£  CljClji^j;^^ 


(2-«) 

PO 

(2-«) 

Po 

(2-«) 

Po 

(2-«) 

Po 

(2-«) 

PO 

(2-«) 

PO 

(2-«) 

PO 

(2-«) 

Po 

(2-«) 

Po 

(2-«) 

PO 

(2-«) 

PO 

(2-«) 

PO 

(2-«) 

Po 

(2-«) 

PO 

(2-«) 

PO 

(2-«) 

PO 

(2-«) 

Po 

^0  °  Lq/D  .  a 

00  °  Lq/D  oo  a 

^0  °  Lq/D  .  a 

^0  °  Lq/D  .  a 

0Q  °  Lo/D  oo  a 

^0  °  Lq/D  .  a 

00  °  Lo/D  .  a 

00  °  Lq/D  oo  a 

00  °  Lo/D  .  a 

00  °  Lo/D  .  a 

0Q  °  Lq/D  oo  a 

00  °  Lo/D  .  a 

00  °  Lo/D  .  a 

00  °  Lo/D  oo  a 

00  °  Lo/D  .  a 

00  °  Lo/D  .  a 

00  °  Lo/D  oo  a 


N 


D, 


ol/ 


VT 


D 


S 


Lo/D  a 


N 

P 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

•  °  <^1!  s 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

■  °  S 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

•  D, 

■  °  S 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

•  °  <^1!  s 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

■  °  S 

N 

p 

Z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

■  °  S 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

•  D, 

°  ^2  / 
^n/  s 

N 

P 

Z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

•  °  <^1!  s 

N 

p 

Z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

■  °  S 

N 

p 

Z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

■  °  S 

N 

P 

z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

z?,  . 

0 

Lo/D 

a 

VT 

•  D, 

°  ^2  / 
^n/  s 

N 

P 

z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

•  °  <^1!  s 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

•  °  S 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

•  °  S 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

Z?,  . 

0 

Lo/D 

a 

VT 

•  D, 

■  °  <^1!  s 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

•  °  S 

N 

p 

Z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

•  °  S 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

Z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

Z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

N 

P 

D,  . 

0 

Lo/D 

a 

VT 

•  D, 

■  °  S 

N 

P 

Z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

N 

P 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

•  °  <^1!  s 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

■  °  S 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

•  D, 

■  °  S 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

•  °  <^1!  s 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

■  °  S 

N 

p 

Z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

■  °  S 

N 

p 

D,  . 

0 

Lo/D 

a 

VT 

•  D, 

°  ^2  / 
^n/  s 

N 

P 

Z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

°  ^2  / 

^n/  s 

N 

P 

D,  . 

0 

Lo/D 

a 

VT 

.  D, 

•  °  <^1!  s 

N 

p 

Z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

■  °  S 

N 

p 

Z?,  . 

0 

Lo/D 

a 

VT 

.  D, 

■  °  S 

N 

P  ■  D, 

0 

Lo/D 

a 

VT 

•  D, 

°  ^2  / 

^n/  s 

N 

P  ■  D, 

0 

Lo/D 

a 

VT 

•  D, 

°  ^2  / 

^n/  s 

N 

P  .  D,  . 

0 

s 

a 

VT 

•  D, 

.  °  • 

N 

P  .  D,  . 

0 

^n/  s 

a 

VT 

•  D, 

•  °  • 

N 

P  .  D,  . 

0 

^2  / 

^n/  s 

a 

VT 

•  D, 

•  °  • 

N 

P  .  D,  . 

0 

^n/  s 

a 

VT 

•  D, 

•  °  • 

N 

P  .  D,  . 

0 

^n/  s 

a 

VT 

•  D, 

•  °  • 

N 

P  .  D,  . 

0 

. 

a 

VT 

•  D, 

•  °  • 

N 

P  .  D,  . 

0 

^2  / 

^n/  s 

a 

VT 

•  D, 

•  °  • 

N 

P  .  D,  . 

0 

^2  / 

^n/  s 

a 

VT 

•  D, 

•  °  Lo/D  . 

N 

P  .  D,  . 

0 

^n/  s 

a 

VT 

•  D, 

•  °  Lo/D  . 

N 

P  .  D,  . 

0 

^n/  s 

a 

VT 

•  D, 

•  °  Lo/D  . 

N 

P  .  D,  . 

0 

. 

a 

VT 

•  D, 

•  °  Lo/D  . 

N 

P  .  D,  . 

0 

^2  / 

^n/  s 

a 

VT 

•  D, 

•  °  Lo/D  . 

N 

P  .  D,  . 

0 

^2  / 

^n/  s 

a 

VT 

• 

•  °  Lo/D 

N 

P 

.  D,  . 

0 

ol/ 

S 

a 

VT 

•  D, 

0 

Lo/D 

N 

P 

.  D,  . 

0 

ol/ 

s 

a 

VT 

•  D, 

0 

Lo/D 

N 

p 

.  D,  . 

0 

^1/ 

s 

a 

VT 

•  D, 

0 

Lo/D 

N 

p 

.  D,  . 

0 

ol/ 

s 

a 

VT 

•  D, 

0 

Lo/D 

N 

p 

.  D,  . 

0 

ol/ 

s 

a 

VT 

•  D, 

0 

Lo/D 

N 

p 

.  D,  . 

0 

ol/ 

s 

a 

VT 

•  D, 

0 

Lo/D 

N 

p 

.  D,  . 

0 

ol/ 

s 

a 

VT 

•  D, 

0 

Lo/D 

N 

p 

.  D,  . 

0 

ol/ 

s 

a 

VT 

•  D, 

0 

Lo/D 

N 

p 

.  D,  . 

0 

ol/ 

s 

a 

VT 

•  D, 

0 

Lo/D 

N 

p 

.  D,  . 

0 

ol/ 

s 

a 

VT 

•  D, 

0 

Lo/D 

N 

p 

.  D,  . 

0 

ol/ 

s 

a 

VT 

•  D, 

0 

Lo/D 

N 

p 

.  D,  . 

0 

ol/ 

s 

a 

VT 

•  D, 

0 

Lo/D  oo 

N 

p 

.  D,  . 

0 

ol/ 

s 

a 

VT 

•  D, 

0 

Lo/D  oo 

N 

p 

.  D,  . 

0 

ol/ 

s 

a 

VT 

•  D, 

0 

Lo/D  oo 

N 

p 

.  D,  . 

0 

ol/ 

s 

a 

VT 

• 

0 

Lo/D  oo 

N 

P  .  D,  . 

0 

^n/  S 

a 

VT 

D,  . 

,  °  Lo/D 

oo 

N 

P  .  D,  . 

0 

^1/  S 

a 

VT 

•  D, 

•  °  Lo/D 

oo 

N 

P 

D, 

0 

VT  JD 

a 

V 

0 

UjD  . 

N 

P 

D, 

0 

VT /D 

a 

V 

0 

Lo/D  . 

N 

P 

D, 

0 

VT /D 

a 

V 

0 

Lo/D  . 

N 

P 

D, 

0 

VT /D 

a 

V 

0 

UID  . 

N 

P 

D, 

0 

VT  JD 

a 

V 

0 

Lo/D  . 

N 

P 

D, 

0 

VT  JD 

a 

V 

0 

Lo/D  . 

N 

P 

D, 

0 

VT /D 

a 

V 

0 

Lo/D  . 

N 

P 

D, 

0 

VT /D 

a 

V 

0 

Lo/D  . 

N 

P 

D, 

0 

VT /D 

a 

V 

0 

Lo/D 

N 

P 

D, 

0 

VT  JD 

a 

V 

0 

Lo/D 

N 

P 

D, 

0 

VT  JD 

a 

V 

0 

Lo/D 

N 

P 

D, 

0 

VT /D 

a 

V 

0 

Lo/D 

N 

P 

D, 

0 

VT /D 

a 

V 

0 

Lo/D 

N 

P 

D, 

0 

vrfD 

a 

V 

°  L^/D 

N 

P 

D, 

0 

vt/D 

a 

V 

°  L^/D 

N 

P 

D, 

0 

vt/D 

a 

V 

°  L^/D 

N 

P 

D, 

0 

vtJD 

a 

V 

°  Lq/D 

oo 

N 

P 

D, 

0 

vrfD 

a 

V 

°  Lq/D 

oo 

N 

p 

D, 

0 

vrfD 

a 

V 

°  UID 

oo 

N 

P 

D, 

0 

vt/D 

a 

V 

°  Lo/D 

oo 

N 

P 

D,  . 

0 

V 

a 

■u 

D 

Lq/D  oo 

N 

P 

D,  . 

0 

V 

a 

■u 

D 

Lq/D  oo 

N 

P 

D,  . 

0 

V 

a 

■u 

D 

Lq/D  oo 

N 

P 

D,  . 

0 

V 

a 

■u 

D 

LqJD  oo 

N 

P 

D,  . 

0 

V 

a 

■u 

D 

LqJD  oo 

N 

P 

D,  . 

0 

V 

a 

■u 

D 

Lq/D  oo 

N 

P 

D,  . 

0 

V 

a 

■u 

D 

Lq/D  oo 

N 

P 

D,  . 

0 

V 

a 

■u 

D 

Lq/D  oo 

N 

P 

•  D, 

0 

V 

a 

■u 

D 

LqJD  oo 

N 

P 

•  D, 

0 

V 

a 

■u 

D 

LqJD  oo 

N 

P 

•  D, 

0 

V 

a 

■u 

D 

Lq/D  oo 

N 

P 

•  D, 

0 

V 

a 

■u 

D 

Lq/D  oo 

N 

P 

•  D, 

0 

V 

a 

■u 

D 

Lq/D  oo 

N 

P 

•  D, 

0 

V 

a 

■u 

D 

LofD  oo 

N 

P 

•  D, 

0 

V 

a 

■u 

D 

LqJD  oo 

N 

P 

•  D, 

0 

V 

a 

■u 

D 

Lq/D  oo 

N 

P 

•  D, 

0 

V 

a 

■u 

D 

Lq/D  oo 

N 

P 

•  D, 

0 

V 

a 

■u 

D 

Lq/D  oo 

N 

P 

• 

0 

V 

a 

■u 

D 

LqJD  oo 

N 

P 

•  D, 

0 

V 

a 

■u 

D 

LqJD  oo 

List  of  Tables 


Zj  ,  ,  f  j  <  rn 


E 

a  jay 

ly 

PO 

2-a 

a 

Lq 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

Lq 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

Lq 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

Lq 

D 

£ 

ajaji 

ly 

Po 

2-a 

a 

Lq 

D 

£ 

ajaji 

ly 

Po 

2-a 

a 

Lq 

D 

£ 

ajaji 

ly 

Po 

2-a 

a 

Lq 

D 

£ 

a  jay 

ly 

Po 

2-a 

a 

Lq 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

Lq 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

Lq 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

Lq 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

Lq 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

Lq 

D 

£ 

ajaji 

ly 

Po 

2-a 

a 

Lq 

D 

£ 

ajaji 

ly 

Po 

2-a 

a 

Lq 

D 

oo 


E 

a  jay 

ly 

Po 

2-a 

a 

11 

3 

^^0 

D 

£ 

a  jay 

ly 

Po 

2-a 

a 

11 

3 

Lq 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

11 

3 

Lq 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

11 

3 

Lq 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

11 

3 

Lq 

D 

oo 

£ 

a  jay 

ly 

PO 

2-a 

a 

Lo 

D 

£ 

ajaji 

ly 

Po 

2-a 

a 

Lo 

D 

£ 

ajaji 

ly 

Po 

2-a 

a 

Lo 

D 

£ 

ajaji 

ly 

Po 

2-a 

a 

Lr 

D 

£ 

a  jay 

ly 

Po 

2-a 

a 

oo 

D  ^ 

£ 

a  jay 

ly 

PO 

2-a 

a 

Lo 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

Lo 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

Lo 

D 

£ 

a  jay 

ly 

PO 

2-a 

a 

Lo 

D 

£ 

ajaji 

ly 

Po 

2-a 

a 

^  OO 

oo 


Abstract 


PERFORMANCE  ANALYSIS 


OF  A 

HARTMANN  WAVEFRONT  SENSOR 
USED  FOR  SENSING 

ATMOSPHERIC  TURBULENCE  STATISTICS 

/.  Introduction 

1.1  Background 

seeing 

ro 

seeing  monitor 


1.2  Atmospheric  Turbulence 


/ 


(j)  h 


(!>  k 


TT 


(3-«) 


a-2 

Po 


Cl 


a 

2 


—  (4  a)  ^2 


2-a 

2 


ko^ 


a 

2 


a-2 


2 


Cl 


a  — 


a  — 


Po 


-  (4-«)7r2A:2a  a  2^  ^  z  z 


a  a  -  (“-4);r-3/2 - 

3-a  ’ 

2 

k  a 


Po 


ro  Cl 


ko  /Lq, 


Lo 


t2  (p  r,  t2 


V  f 

ti  <t2  (p  f,  ti 


(f)  r,t2  (f)  r  V  t2  -  h  ,ti  . 


P^i 


Vi 


Q 


P^tot 


p^i 


2-a 

POtot 


Q 


2-a 
P^i  ■ 


1  =  1 


1.3  Hartman  Wavefront  Sensor  as  Turbulence  Monitor 


Ds  p,qMH2 


£  Sa  P,  h  -  Sj  q,  t2 


t2 


p 


Sa  P,  h 


h 

a  Sj  q,t2 


q 


b 


p,q,ti,t2  £  sjpHi  sf  q,t2  -  Sa  pM  s^  qM 


5 


s,^q-p  vt2-h 


£  Sa  p,  h  Sj  q, t2 


a  b 


tion  a  b 

function 


self  slope  structure  func- 
cross  slope  structure 


a 


b 


sQ-p  vt2-ti  £  {s  p,ti  s  q,t2  }  . 


p  p-q 


P 


M 


p  MN 

P 


M  N 


Ds  p 


NM 


•S  Qmj  tn  S  Pm:  t-n 


m=l  n=l 


Ds  p 


P  Qm  Pm 
Pm 


P  qm  -  Pm  rn 
•S  Pm:  tfi 


s 


SNR  Ds  p 


£  Ds  p  - 


2  2  ’ 
£  Ds  p  -  £  Ds  p 


Ds 


1-4  Problem  Description  and  Scope 


1.5  Organization 


II.  Slope  Correlation  Funetion  and  Slope  Strueture  Funetion  (SSF) 


S 


2.1  Zernike  Polynomials 


3 


Z 

Z 


r 

■‘‘■n  ' 


3 

r  y/ 

m9 

m  / 

3 

R^  r 

mO 

Zj  r,  e 

~Kr 

m 

(n—m)/2 

—  ^  n  —  s 

„n-2s 

.=0  ^ 

n 

ml  —  s  n  — 

■  m  / 

f 

—  s 

n  m 


n 


m  <n 
n  —  \m 


m 


00 


r  6W  r  Zj  r,  6  Zjt  r,  9  5jji , 


R 


W  r 


/tt  r  < 


r  > 


3  3 


(f>  r,  6 


cj)  Rp,6  0‘j^j  5 

j 


P  r/R, 


LXJ 

“i  r  9W  -  (f>  r,e  Zj  -,9 

JL  — cso  -tX,  -tX, 


9^  j  k  ^  (j) 


Zj  p,9 


W  p  Zj  p,  9 


Qj  k,(f)  —  irik  ■  p  , 


-  rn4> 

Q  j  k,4> 


2  ^  5  ^ 


J 

j  ^  5  ^ 


TT 


_  n/2 


m 


2.2  Zernike  Expansion  Coefficient  Spatial  Covariance  for  Turbulence  Induced  Phase 


D  DR 

xy 

u  u  u 

6o  u  X 

00 

(f>  Rp  P  5 

i=i 

oo 


(f>  R  p'  u  Oji  u  Zji  p'  . 


00 

Oj  (f)  Rp  Zj  p  W  p  p, 

—  00 


00 

Oji  u  (f)  R  p'  u  Zji  p'  W  p'  p' , 

—  00 

W  p 


00  00 

£  OjO*!  u  p  p'£  (f)  Rp  (jf  R  p'  u  Zj  p  Z*i  p'  W  p  W  p'  , 

—  00  —00 


oo  27r 

£  a,a*,  u  k  e  -i2n2ukcos{0-0o)  ,  k/R  Qi  k,e  Q*,  k,e  , 


Qj  k,0 


k 


Po 


ciz 


£  ajG*/  u 


Cl 


]J  («-2) 

Po 


a 

2 


2-a 

2 


n  n 


1/2 


fij'  W,  00,^0 


fjji  u,0o,ko 


m,n' /  j,j' 

-  (r^+n'-m+m')/2  ^  m'  6^  Im+m' ,n+l,n'+l  U,ko 

_  ^n+n'+2m+\m-m'\)/2  m  -  m'  0o  /|„' |,n+l,n'+l  ^’^0 

m,n' /  j,j' 

__  ^n+n'-m+m')/2  ^  m'  Oo  Im+m' ,n+l,n'+l  U,ko 

_  (n+n'+2m+|m-m'|)/2  m  —  m'  ftr.  T,  ,,  ,  ,  ,  11  i-n 

m,  n'  /  j  j' 

-  '(-+-'-™+™')/2  ^  m'  00  Wm',n+l,n'+l  ^0 

__  in+n'+2m+\m-m'\)l2  m  -  m'  0o  /|„' |,n+l,n'+l  ^’^0 

m,  n'  /  j  j' 

-  'in+n'-m+m')/2  rn  m'  6^  n  Im+m',n+l,n'+l  U,ko 

__  (n+n'+2m+\m-m'\)/2  rn  -  m!  0^  TT  uM 


m  J 


- 

{n+n' -m)/2  ^ 

1X100  Im,n+l,n'+l  ^0 

m 

f 

(n+n' -m)/2^ 

— 

m  00 

1^  -^m,n+l,n'+l 

u,ko 

m' 

j 

{n+n' -m)/2  ^ 

— 

1X100  Im,n+l,n'+l  H:  ^0 

m 

f 

{n+n' -m)/2  ^ 

— 

m  00 

1^  -^m,n+l,n'+l 

u,ko 

m 

m' 

("+"')/27o,n+l,n'+l  ^>^0 

D 

1  0 

0 

J 

a,  xo 

°°  X  Vk  ax  Jn 

X  Jp 

0 

to 

0 

to 

a/2 

X, 


2.3  Slope  Correlation  Function  and  Slope  Structure  Function  Evaluation 


S 


Ds 


Sa  X,t 


rW  r  —  X  Wcj)  r,  t 


Sf,  X,  t 


a 


t 

X  W  r 

(p  r,t  V 

W  r 


s  x,t 


X 


et.  al. 


^Evaluating  Eqn.  (2.21)  is  computational  expensive  and  ultimately  limits  which  problems  can  realisti¬ 
cally  be  solved.  Reference  [15]  presents  a  closed  form  solution  for  Eqn.  (2.21),  where  u  >  2,  consisting  of 
infinite  sums  of  hypergeometric  functions  of  type  4F3.  They  also  state  it  is  easier  to  numerically  integrate 
Eqn.  (2.21)  than  to  use  the  closed  form.  This  author  has  explored  both  approaches  and  agrees. 
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V.  Results,  Conclusions  and  Recommendations 
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5.3  Recommendation  for  Future  Research 


^This  simulation  should  include  a  phase  screen  generator  capable  of  including  outer  scale,  power  law,  and 
temporal  effects.  This  author  has  developed  and  tested  a  phase  screen  generator,  with  all  these  capabilities, 
based  on  Eqn.  (2.18). 


Appendix  A.  Zernike  Expansion  Coefficient  Covariance  Plots 


S 


S  Gja*/  u 


]J  («-2) 
Cl  — 

Po 


n 


n 


1/2 


fjji  u,Oo,ko 


fjj’  u,6o,ko 

_  (n+n')/2  00  l2,n+l,n'+l  uffio  j,j' 

-  n,ko 

__  in+n')/2  00  72,n+l,n'+l  ^0  J,  / 

-  n,ko 


Cl 


00 


ko  71  —  , 

-^0 


O’-!  ^0 


X  ax  Jfj  X  Ju  X 

9  9  q;/2 

X‘^  Xo 


X, 


D 


Ln 


J X  Jn  X  J  1/  X 


K  fj,  V 


U 


a 


LoD 


Appendix  B.  Zernike  Expansion  Coefficient  Covariance  Tables 
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Appendix  C.  SSF  Estimator  Second  Moment  Simplification 


(* - 

This  Mathematica  code  helps  determine  the  second  moment  for  the 
slope  structure  function  estimator. 


x2-xl=p 

x4-x3=p 


easel:  xl=x3  and  x2=x4  and  delta_time=0 
case2:  x2=x3  and  delta_time=0 
case3:  x4=xl  and  delta_time=0 

case4:  xl,x2,x3,and  x4  are  four  unique  location  (no  overlap) 

Toby  D.  Reeves,  July-Sept  1996 
*) 

- 

*) 

Remove [ " Global ‘ ] 

(* - 

Define  the  Expected  Value  Operator 
*) 

EV [a_+b_] : =EV [a] +EV [b] 

EV[a_Integer  b_] :=a  EV[b] 

- 

Define  a  Wrapper  zmgrv[  ]  for  Zero  Mean  Gaussian  Random  Variable, 
then  define  rules  that  apply  to  zmgrv. 

*) 

Format [zmgrv [x_] ] : =x 

EV[ul_ zmgrv  u2_ zmgrv  u3_ zmgrv  u4_ zmgrv] :=  ( 

EVEul  u2]  EV[u3  u4]  + 

EVEul  u3]  EVEu2  u4]  + 

EVEul  u4]  EVEu2  u3] 

) 

EVEul_ zmgrv  ul_zmgrv  u3_ zmgrv  u4_ zmgrv] :=  ( 

EVEul  u2]  EVEu3  u4]  + 


EV[ul  u3]  EV[u2  u4]  + 

EV[ul  u4]  EV[u2  u3] 

) / .u2->ul 

EV[ul_zmgrv  ul_zmgrv  u3_zmgrv  u3_zmgrv] :=  ( 

EV[ul  u2]  EV[u3  u4]  + 

EV[ul  u3]  EV[u2  u4]  + 

EV[ul  u4]  EV[u2  u3] 

) / . {u2->ul ,u4->u3} 

EV[ul_zmgrv  ul_zmgrv  ul_zmgrv  u4_zmgrv] :=  ( 

EV[ul  u2]  EV[u3  u4]  + 

EV[ul  u3]  EV[u2  u4]  + 

EV[ul  u4]  EV[u2  u3] 

) / . {u2->ul ,u3->ul} 

EV[ul_zmgrv  ul_zmgrv  ul_zmgrv  ul_zmgrv] :=  ( 

EV[ul  u2]  EV[u3  u4]  + 

EV[ul  u3]  EV[u2  u4]  + 

EV[ul  u4]  EV[u2  u3] 

) / . {u2->ul ,u3->ul ,u4->ul} 

- 

Properties  of  the  Noise 
*) 

noiseRules=  { 

EV [zmgrv [n [a_]  ]  zmgrv [n [a_]  ]  ] ->sig"2 , 

EV [zmgrv [n [a_]  ]  zmgrv [n [b_] ] ] ->0 , 

EV [zmgrv [n [a_] ] ] ->0 , 

EV [zmgrv [n [a_] ]  zmgrv [s [x_] ]] ->0 

>; 

(* - 

Apply  these  rules  to  my  problem  of  calculating  the  expectation 
within  the  Second  Moment  formula 

*) 

fl=  (a-b)"2  (c-d)"2  /. 

{a->zmgrv [s [xl] ] +zmgrv [n [xl] ] , 
b->zmgrv [s [x2] ] +zmgrv [n [x2] ] , 
c->zmgrv [s [x3] ] +zmgrv [n [x3] ] , 
d->zmgrv [s [x4] ] +zmgrv [n [x4]  ] 

>; 

BigMess=EV [Expand [fl]]  ; 
convertRule={ 

EV[zmgrv[s[x_]]  "2]  ->G[0], 

EV [zmgrv [s[xl]]  zmgrv [s [x2] ] ]  ->G[x2-xl], 

EV [zmgrv [s  [x3] ]  zmgrv [s  [x4] ] ]  ->G [x4-x3] , 


EV [zmgrv [s[xl]]  zmgrv [s [x3] ] ] 


->G  [x3-xl] , 


EV[zmgrv[s[x2]]  zmgrv  [s  [x4]  ]  ] 


->G [x4-x2]  , 


EV [zmgrv  [s  [xl] ]  zmgrv [s  [x4] ] ]  ->G[x4-xl], 

EV [zmgrv  [s  [x2] ]  zmgrv [s  [x3] ] ]  ->G[x3-x2] 


(*  apply  geometry  and  noiseRules  for  each  case  *) 

caselRule={xl->x3,x2->x4} ; 

casela= (BigMess/ . caselRule) / .noiseRules ; 

case2Rule={x2->x3} ; 

case2a= (BigMess/ . case2Rule) / .noiseRules ; 
case3Rule={xl->x4} ; 

case3a= (BigMess/ . case3Rule) / . noiseRules ; 
case4Rule={} ; 

case4a= (BigMess/ . case4Rule) / . noiseRules ; 

(*  write  in  Gamma  notation  *) 
caselb=casela/ . convertRule ; 
case2b=case2a/ . convertRule ; 
case3b=case3a/ . convertRule ; 
case4b=case4a/ . convertRule ; 

(*  use  geometry  to  simplify  more  *) 

geoRule={x2-xl->p,x4-x3->p,x4-x2->x3-xl}  (♦valid  all  cases*) 
caselc=ExpandAll [caselb/ . geoRule] ; 

case2c=ExpandAll [(case2b/ . geoRule) / . {x3-xl->p,x4-xl->2p}]  ; 
case3c=ExpandAll [(case3b/ . geoRule) / . {x3-xl->-p,x3-x2->-2p}] ; 
case4c=ExpandAll [case4b/ . geoRule] ; 

(*  We  know  that  for  our  case  G[p]=G[-p]  *) 
q=case3c 

case3c=case3c/ . {-2p  ->  2p,  -p  ->  p}; 

- 

Drum  roll  please! I 
*) 

sout=OpenWrite ["sf5n.txt" ,FormatType->0utputForm,PageWidth->75] 

Write [sout , "BigMess=  " .ExpandAll [BigMess] ] 

Write [sout , "casela=  ",casela] 

Write [sout , "case2a=  ",case2a] 


Write [sout , "case3a=  ",case3a] 

Write [sout , "case4a=  ",case4a] 

Write [sout , "caselb=  ",caselb] 

Write [sout , "case2b=  ",case2b] 

Write [sout , "case3b=  ",case3b] 

Write [sout , "case4b=  ",case4b] 

Write [sout , "caselc=  ",caselc] 

Write [sout , "case2c=  ",case2c] 

Write [sout , "case3c=  ",case3c] 

Write [sout , "case4c=  ",case4c] 

Close [sout] 

2 

BigMess=  2  EV[n[xl]  n[x3]]  -  4  EV[n[xl]  n[x3]]  EV[n[x2]  n[x3]]  + 

2  2  2 

>  2  EV[n[x2]  n[x3]]  +  EV[n[xl]  ]  EV[n[x3]  ]  - 

2  2  2 

>  2  EV[n[xl]  n[x2]]  EV[n[x3]  ]  +  EV[n[x2]  ]  EV[n[x3]  ]  - 

>  4  EV[n[xl]  n[x3]]  EV[n[xl]  n[x4]]  + 

2 

>  4  EV[n[x2]  n[x3]]  EV[n[xl]  n[x4]]  +  2  EV[n[xl]  n[x4]]  + 

>  4  EV[n[xl]  n[x3]]  EV[n[x2]  n[x4]]  - 

>  4  EV[n[x2]  n[x3]]  EV[n[x2]  n[x4]]  - 

2 

>  4  EV[n[xl]  n[x4]]  EV[n[x2]  n[x4]]  +  2  EV[n[x2]  n[x4]] 

2 

>  2  EV[n[xl]  ]  EV[n[x3]  n[x4]]  +  4  EV[n[xl]  n[x2]]  EV[n[x3]  n[x4]]  - 

2  2  2 

>  2  EV[n[x2]  ]  EV[n[x3]  n[x4]]  +  EV[n[xl]  ]  EV[n[x4]  ]  - 

2  2  2 

>  2  EV[n[xl]  n[x2]]  EV[n[x4]  ]  +  EV[n[x2]  ]  EV[n[x4]  ]  + 


2 


>  2  EV[n[x3]  ]  EV[n[xl]  s[xl]]  -  4  EV[n[x3]  n[x4]]  EV[n[xl]  s[xl]]  + 

2  2 

>  2  EV[n[x4]  ]  EV[n[xl]  s[xl]]  -  2  EV[n[x3]  ]  EV[n[x2]  s[xl]]  + 

2 

>  4  EV[n[x3]  n[x4]]  EV[n[x2]  s[xl]]  -  2  EV[n[x4]  ]  EV[n[x2]  s[xl]]  + 

>  4  EV[n[xl]  n[x3]]  EV[n[x3]  s[xl]]  - 

>  4  EV[n[x2]  n[x3]]  EV[n[x3]  s[xl]]  - 

>  4  EV[n[xl]  n[x4]]  EV[n[x3]  s[xl]]  + 

2 

>  4  EV[n[x2]  n[x4]]  EV[n[x3]  s[xl]]  +  2  EV[n[x3]  s[xl]] 

>  4  EV[n[xl]  n[x3]]  EV[n[x4]  s[xl]]  + 

>  4  EV[n[x2]  n[x3]]  EV[n[x4]  s[xl]]  + 

>  4  EV[n[xl]  n[x4]]  EV[n[x4]  s[xl]]  - 

>  4  EV[n[x2]  n[x4]]  EV[n[x4]  s[xl]]  - 

2 

>  4  EV[n[x3]  s[xl]]  EV[n[x4]  s[xl]]  +  2  EV[n[x4]  s[xl]]  + 

2  2  2 

>  EV[n[x3]  ]  EV[s[xl]  ]  -  2  EV[n[x3]  n[x4]]  EV[s[xl]  ]  + 

2  2  2 

>  EV[n[x4]  ]  EV[s[xl]  ]  -  2  EV[n[x3]  ]  EV[n[xl]  s[x2]]  + 

2 

>  4  EV[n[x3]  n[x4]]  EV[n[xl]  s[x2]]  -  2  EV[n[x4]  ]  EV[n[xl]  s[x2]]  + 

2 

>  2  EV[n[x3]  ]  EV[n[x2]  s[x2]]  -  4  EV[n[x3]  n[x4]]  EV[n[x2]  s[x2]]  + 

2 

>  2  EV[n[x4]  ]  EV[n[x2]  s[x2]]  -  4  EV[n[xl]  n[x3]]  EV[n[x3]  s[x2]]  + 

>  4  EV[n[x2]  n[x3]]  EV[n[x3]  s[x2]]  + 


> 

4 

EV  [n[xl] 

n[x4]] 

EV [n [x3] 

s[x2]]  - 

> 

4 

EV  [n[x2] 

n[x4]] 

EV  [n[x3] 

s[x2]]  - 

> 

4 

EV [n [x3] 

s[xl]] 

EV  [n[x3] 

s [x2] ]  + 

> 

4 

EV  [n[x4] 

s[xl]] 

EV  [n[x3] 

2 

s[x2]]  +  2  EV[n[x3]  s[x2]]  + 

> 

4 

EV  [n[xl] 

n[x3]] 

EV  [n[x4] 

s[x2]]  - 

> 

4 

EV  [n[x2] 

n[x3]] 

EV  [n[x4] 

s[x2]]  - 

> 

4 

EV  [n[xl] 

n[x4]] 

EV  [n[x4] 

s [x2] ]  + 

> 

4 

EV  [n[x2] 

n[x4]] 

EV  [n[x4] 

s [x2] ]  + 

> 

4 

EV  [n[x3] 

s[xl]] 

EV  [n[x4] 

s[x2]]  - 

> 

4 

EV  [n[x4] 

s[xl]] 

EV  [n[x4] 

s[x2]]  - 

> 

4 

EV [n [x3] 

s[x2]] 

EV  [n[x4] 

2 

s[x2]]  +  2  EV[n[x4]  s  [x2]  ] 

> 

2 

EV  [n[x3] 

2 

]  EV[s[xl]  s[x2]]  +  4  EV[n[x3]  n[x4]]  EV[s[xl]  s[x2]]  - 

> 

2 

EV  [n[x4] 

2  2  2 

]  EV[s[xl]  s[x2]]  +  EV[n[x3]  ]  EV[s[x2]  ]  - 

> 

2 

EV [n [x3] 

n[x4]] 

EV[s[x2] 

2  2  2 

]  +  EV[n[x4]  ]  EV[s[x2]  ]  + 

> 

4 

EV  [n[xl] 

n[x3]] 

EV  [n[xl] 

s [x3] ]  - 

> 

4 

EV  [n[x2] 

n[x3]] 

EV  [n[xl] 

s [x3] ]  - 

> 

4 

EV  [n[xl] 

n[x4]] 

EV  [n[xl] 

s [x3] ]  + 

> 

4 

EV  [n[x2] 

n[x4]] 

EV  [n[xl] 

s [x3] ]  + 

> 

4 

EV  [n[x3] 

s[xl]] 

EV  [n[xl] 

s [x3] ]  - 

> 

4 

EV  [n[x4] 

s[xl]] 

EV  [n[xl] 

s [x3] ]  - 

> 


4  EV[n[x3]  s[x2]]  EV[n[xl]  s[x3]]  + 


2 

>  4  EV[n[x4]  s[x2]]  EV[n[xl]  s[x3]]  +  2  EV[n[xl]  s[x3]] 

>  4  EV[n[xl]  n[x3]]  EV[n[x2]  s[x3]]  + 

>  4  EV[n[x2]  n[x3]]  EV[n[x2]  s[x3]]  + 

>  4  EV[n[xl]  n[x4]]  EV[n[x2]  s[x3]]  - 

>  4  EV[n[x2]  n[x4]]  EV[n[x2]  s[x3]]  - 

>  4  EV[n[x3]  s[xl]]  EV[n[x2]  s[x3]]  + 

>  4  EV[n[x4]  s[xl]]  EV[n[x2]  s[x3]]  + 

>  4  EV[n[x3]  s[x2]]  EV[n[x2]  s[x3]]  - 

>  4  EV[n[x4]  s[x2]]  EV[n[x2]  s[x3]]  - 

2 

>  4  EV[n[xl]  s[x3]]  EV[n[x2]  s[x3]]  +  2  EV[n[x2]  s[x3]]  + 

2 

>  2  EV[n[xl]  ]  EV[n[x3]  s[x3]]  -  4  EV[n[xl]  n[x2]]  EV[n[x3]  s[x3]]  + 

2 

>  2  EV[n[x2]  ]  EV[n[x3]  s[x3]]  +  4  EV[n[xl]  s[xl]]  EV[n[x3]  s[x3]]  - 

2 

>  4  EV[n[x2]  s[xl]]  EV[n[x3]  s[x3]]  +  2  EV[s[xl]  ]  EV[n[x3]  s[x3]]  - 

>  4  EV[n[xl]  s[x2]]  EV[n[x3]  s[x3]]  + 

>  4  EV[n[x2]  s[x2]]  EV[n[x3]  s[x3]]  - 

2 

>  4  EV[s[xl]  s[x2]]  EV[n[x3]  s[x3]]  +  2  EV[s[x2]  ]  EV[n[x3]  s[x3]]  - 

2 

>  2  EV[n[xl]  ]  EV[n[x4]  s[x3]]  +  4  EV[n[xl]  n[x2]]  EV[n[x4]  s[x3]]  - 

2 

>  2  EV[n[x2]  ]  EV[n[x4]  s[x3]]  -  4  EV[n[xl]  s[xl]]  EV[n[x4]  s[x3]]  + 


2 

>  4  EV[n[x2]  s[xl]]  EV[n[x4]  s[x3]]  -  2  EV[s[xl]  ]  EV[n[x4]  s[x3]]  + 

>  4  EV[n[xl]  s[x2]]  EV[n[x4]  s[x3]]  - 

>  4  EV[n[x2]  s[x2]]  EV[n[x4]  s[x3]]  + 

2 

>  4  EV[s[xl]  s[x2]]  EV[n[x4]  s[x3]]  -  2  EV[s[x2]  ]  EV[n[x4]  s[x3]]  + 

>  4  EV[n[xl]  n[x3]]  EV[s[xl]  s[x3]]  - 

>  4  EV[n[x2]  n[x3]]  EV[s[xl]  s[x3]]  - 

>  4  EV[n[xl]  n[x4]]  EV[s[xl]  s[x3]]  + 

>  4  EV[n[x2]  n[x4]]  EV[s[xl]  s[x3]]  + 

>  4  EV[n[x3]  s[xl]]  EV[s[xl]  s[x3]]  - 

>  4  EV[n[x4]  s[xl]]  EV[s[xl]  s[x3]]  - 

>  4  EV[n[x3]  s[x2]]  EV[s[xl]  s[x3]]  + 

>  4  EV[n[x4]  s[x2]]  EV[s[xl]  s[x3]]  + 

>  4  EV[n[xl]  s[x3]]  EV[s[xl]  s[x3]]  - 

2 

>  4  EV[n[x2]  s[x3]]  EV[s[xl]  s[x3]]  +  2  EV[s[xl]  s[x3]] 

>  4  EV[n[xl]  n[x3]]  EV[s[x2]  s[x3]]  + 

>  4  EV[n[x2]  n[x3]]  EV[s[x2]  s[x3]]  + 

>  4  EV[n[xl]  n[x4]]  EV[s[x2]  s[x3]]  - 

>  4  EV[n[x2]  n[x4]]  EV[s[x2]  s[x3]]  - 

>  4  EV[n[x3]  s[xl]]  EV[s[x2]  s[x3]]  + 

>  4  EV[n[x4]  s[xl]]  EV[s[x2]  s[x3]]  + 


> 


4  EV[n[x3]  s[x2]]  EV[s[x2]  s[x3]] 


> 


4  EV[n[x4]  s[x2]]  EV[s[x2]  s  [x3]  ] 


>  4  EV[n[xl]  s[x3]]  EV[s[x2]  s[x3]]  + 

>  4  EV[n[x2]  s[x3]]  EV[s[x2]  s[x3]]  - 

2 

>  4  EV[s[xl]  s[x3]]  EV[s[x2]  s[x3]]  +  2  EV[s[x2]  s[x3]]  + 

2  2  2 

>  EV[n[xl]  ]  EV[s[x3]  ]  -  2  EV[n[xl]  n[x2]]  EV[s[x3]  ]  + 

2  2  2 

>  EV[n[x2]  ]  EV[s[x3]  ]  +  2  EV[n[xl]  s[xl]]  EV[s[x3]  ]  - 

2  2  2 

>  2  EV[n[x2]  s[xl]]  EV[s[x3]  ]  +  EV[s[xl]  ]  EV[s[x3]  ]  - 

2  2 

>  2  EV[n[xl]  s[x2]]  EV[s[x3]  ]  +  2  EV[n[x2]  s[x2]]  EV[s[x3]  ]  - 

2  2  2 

>  2  EV[s[xl]  s[x2]]  EV[s[x3]  ]  +  EV[s[x2]  ]  EV[s[x3]  ]  - 

>  4  EV[n[xl]  n[x3]]  EV[n[xl]  s[x4]]  + 

>  4  EV[n[x2]  n[x3]]  EV[n[xl]  s[x4]]  + 

>  4  EV[n[xl]  n[x4]]  EV[n[xl]  s[x4]]  - 

>  4  EV[n[x2]  n[x4]]  EV[n[xl]  s[x4]]  - 

>  4  EV[n[x3]  s[xl]]  EV[n[xl]  s[x4]]  + 

>  4  EV[n[x4]  s[xl]]  EV[n[xl]  s[x4]]  + 

>  4  EV[n[x3]  s[x2]]  EV[n[xl]  s[x4]]  - 

>  4  EV[n[x4]  s[x2]]  EV[n[xl]  s[x4]]  - 

>  4  EV[n[xl]  s[x3]]  EV[n[xl]  s[x4]]  + 


> 


4  EV[n[x2]  s[x3]]  EV[n[xl]  s[x4]] 


> 

4 

EV[s[xl] 

s [x3] ] 

EV  [n[xl] 

s  [x4]  ] 

+ 

> 

4 

EV[s[x2] 

s [x3] ] 

EV  [n[xl] 

s  [x4]  ] 

2 

+  2  EV[n[xl]  s[x4]]  + 

> 

4 

EV  [n[xl] 

n[x3]] 

EV  [n[x2] 

s  [x4]  ] 

- 

> 

4 

EV  [n[x2] 

n[x3]] 

EV  [n[x2] 

s  [x4]  ] 

- 

> 

4 

EV  [n[xl] 

n[x4]] 

EV  [n[x2] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[x2] 

n[x4]] 

EV  [n[x2] 

s  [x4]  ] 

+ 

> 

4 

EV [n [x3] 

s[xl]] 

EV  [n[x2] 

s  [x4]  ] 

- 

> 

4 

EV  [n[x4] 

s[xl]] 

EV  [n[x2] 

s  [x4]  ] 

- 

> 

4 

EV  [n[x3] 

s[x2]] 

EV  [n[x2] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[x4] 

s[x2]] 

EV  [n[x2] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[xl] 

s [x3] ] 

EV  [n[x2] 

s  [x4]  ] 

- 

> 

4 

EV  [n[x2] 

s [x3] ] 

EV  [n[x2] 

s  [x4]  ] 

+ 

> 

4 

EV[s[xl] 

s [x3] ] 

EV  [n[x2] 

s  [x4]  ] 

- 

> 

4 

EV[s[x2] 

s [x3] ] 

EV  [n[x2] 

s  [x4]  ] 

- 

> 

4 

EV  [n[xl] 

s [x4] ] 

EV  [n[x2] 

s  [x4]  ] 

2 

+  2  EV[n[x2]  s[x4]] 

> 

2 

EV  [n[xl] 

2 

]  EV[n[x3]  s[x4]]  +  4 

EV[n[xl]  n[x2]]  EV[n[x3] 

s  [x4]  ]  - 

> 

2 

EV  [n[x2] 

2 

]  EV[n[x3]  s[x4]]  -  4 

EV[n[xl]  s[xl]]  EV[n[x3] 

s  [x4]  ]  + 

> 

4 

EV  [n[x2] 

s[xl]] 

EV  [n[x3] 

s  [x4]  ] 

2 

-  2  EV[s[xl]  ]  EV[n[x3] 

s  [x4]  ]  + 

> 

4 

EV  [n[xl] 

s[x2]] 

EV  [n[x3] 

s  [x4]  ] 

- 

> 

4 

EV  [n[x2] 

s[x2]] 

EV [n [x3] 

s  [x4]  ] 

+ 

2 


> 

4 

EV[s[xl] 

s[x2]] 

EV  [n[x3] 

s  [x4]  ] 

-  2  EV[s[x2]  ] 

EV  [n[x3] 

s  [x4]  ]  + 

> 

2 

EV  [n[xl] 

2 

]  EV[n[x4]  s[x4]]  -  4 

EV[n[xl]  n[x2]] 

EV  [n[x4] 

s  [x4]  ]  + 

> 

2 

EV  [n[x2] 

2 

]  EV[n[x4]  s[x4]]  +  4 

EV[n[xl]  s[xl]] 

EV  [n[x4] 

s  [x4]  ]  - 

> 

4 

EV  [n[x2] 

s[xl]] 

EV  [n[x4] 

s  [x4]  ] 

2 

+  2  EV[s[xl]  ] 

EV  [n[x4] 

s  [x4]  ]  - 

> 

4 

EV  [n[xl] 

s[x2]] 

EV  [n[x4] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[x2] 

s[x2]] 

EV  [n[x4] 

s  [x4]  ] 

- 

> 

4 

EV[s[xl] 

s[x2]] 

EV  [n[x4] 

s  [x4]  ] 

2 

+  2  EV[s[x2]  ] 

EV  [n[x4] 

s  [x4]  ]  - 

> 

4 

EV  [n[xl] 

n[x3]] 

EV[s[xl] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[x2] 

n[x3]] 

EV[s  [xl] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[xl] 

n[x4]] 

EV[s  [xl] 

s  [x4]  ] 

- 

> 

4 

EV  [n[x2] 

n[x4]] 

EV[s[xl] 

s  [x4]  ] 

- 

> 

4 

EV  [n[x3] 

s[xl]] 

EV[s[xl] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[x4] 

s[xl]] 

EV[s[xl] 

s  [x4]  ] 

+ 

> 

4 

EV [n [x3] 

s[x2]] 

EV[s  [xl] 

s  [x4]  ] 

- 

> 

4 

EV  [n[x4] 

s[x2]] 

EV[s[xl] 

s  [x4]  ] 

- 

> 

4 

EV  [n[xl] 

s [x3] ] 

EV[s  [xl] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[x2] 

s [x3] ] 

EV[s[xl] 

s  [x4]  ] 

- 

> 

4 

EV[s  [xl] 

s [x3] ] 

EV[s[xl] 

s  [x4]  ] 

+ 

> 

4 

EV[s[x2] 

s [x3] ] 

EV[s[xl] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[xl] 

s [x4] ] 

EV[s[xl] 

s  [x4]  ] 

_ 

> 

4 

EV[n[x2] 

s [x4] ] 

EV[s[xl] 

s  [x4]  ] 

2 

+  2  EV[s[xl]  s[x4]]  + 

> 

4 

EV  [n[xl] 

n[x3]] 

EV[s[x2] 

s  [x4]  ] 

- 

> 

4 

EV  [n[x2] 

n[x3]] 

EV[s[x2] 

s  [x4]  ] 

- 

> 

4 

EV  [n[xl] 

n[x4]] 

EV[s[x2] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[x2] 

n[x4]] 

EV[s[x2] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[x3] 

s[xl]] 

EV[s[x2] 

s  [x4]  ] 

- 

> 

4 

EV  [n[x4] 

s[xl]] 

EV[s[x2] 

s  [x4]  ] 

- 

> 

4 

EV  [n[x3] 

s[x2]] 

EV[s[x2] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[x4] 

s[x2]] 

EV[s[x2] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[xl] 

s [x3] ] 

EV[s[x2] 

s  [x4]  ] 

- 

> 

4 

EV  [n[x2] 

s [x3] ] 

EV[s[x2] 

s  [x4]  ] 

+ 

> 

4 

EV[s  [xl] 

s [x3] ] 

EV[s[x2] 

s  [x4]  ] 

- 

> 

4 

EV[s[x2] 

s [x3] ] 

EV[s[x2] 

s  [x4]  ] 

- 

> 

4 

EV  [n[xl] 

s [x4] ] 

EV[s[x2] 

s  [x4]  ] 

+ 

> 

4 

EV  [n[x2] 

s [x4] ] 

EV[s[x2] 

s  [x4]  ] 

- 

> 

4 

EV[s[xl] 

s [x4] ] 

EV[s[x2] 

s  [x4]  ] 

2 

+  2  EV[s[x2]  s[x4]] 

> 

2 

EV  [n[xl] 

2 

]  EV[s[x3]  s[x4]]  +  4 

EV[n[xl]  n[x2]]  EV[s[x3] 

s  [x4]  ]  - 

> 

2 

EV  [n[x2] 

2 

]  EV[s[x3]  s[x4]]  -  4 

EV[n[xl]  s[xl]]  EV[s[x3] 

s  [x4]  ]  + 

> 

4 

EV  [n[x2] 

s[xl]] 

EV[s[x3] 

s  [x4]  ] 

2 

-  2  EV[s[xl]  ]  EV[s[x3] 

s  [x4]  ]  + 

> 

4 

EV  [n[xl] 

s[x2]] 

EV[s[x3] 

s  [x4]  ] 

_ 

>  4  EV[n[x2]  s[x2]]  EV[s[x3]  s[x4]]  + 


2 

>  4  EV[s[xl]  s[x2]]  EV[s[x3]  s[x4]]  -  2  EV[s[x2]  ]  EV[s[x3]  s[x4]]  + 

2  2  2 

>  EV[n[xl]  ]  EV[s[x4]  ]  -  2  EV[n[xl]  n[x2]]  EV[s[x4]  ]  + 

2  2  2 

>  EV[n[x2]  ]  EV[s[x4]  ]  +  2  EV[n[xl]  s[xl]]  EV[s[x4]  ]  - 

2  2  2 

>  2  EV[n[x2]  s[xl]]  EV[s[x4]  ]  +  EV[s[xl]  ]  EV[s[x4]  ]  - 

2  2 

>  2  EV[n[xl]  s[x2]]  EV[s[x4]  ]  +  2  EV[n[x2]  s  [x2]  ]  EV[s[x4]  ]  - 

2  2  2 

>  2  EV[s[xl]  s[x2]]  EV[s[x4]  ]  +  EV[s[x2]  ]  EV[s[x4]  ] 

4  2  2  2  2 

casela=  12  sig  +  12  sig  EV[s[x3]  ]  +  3  EV[s[x3]  ] 

2  2 

>  24  sig  EV[s[x3]  s  [x4]  ]  -  12  EV[s[x3]  ]  EV[s[x3]  s[x4]]  + 

2  2  2  2  2 

>  4  EV[s[x3]  s[x4]]  +  12  sig  EV[s[x4]  ]  +  2  EV[s[x3]  ]  EV[s[x4]  ]  - 

2  2  2 

>  12  EV[s[x3]  s[x4]]  EV[s[x4]  ]  +  3  EV[s[x4]  ]  + 

2  2  2 

>  4  (2  EV[s[x3]  s[x4]]  +  EV[s[x3]  ]  EV[s[x4]  ]) 

4  2  2  2 

case2a=  6  sig  +  2  sig  EV[s[xl]  ]  -  8  sig  EV[s[xl]  s[x3]]  + 

2  2  2  2  2 

>  2  EV[s[xl]  s[x3]]  +  8  sig  EV[s[x3]  ]  +  EV[s[xl]  ]  EV[s[x3]  ]  - 

2  2  2 

>  6  EV[s[xl]  s[x3]]  EV[s[x3]  ]  +  3  EV[s[x3]  ]  + 

2  2 

>  4  sig  EV[s[xl]  s[x4]]  +  2  EV[s[xl]  s  [x4]  ] 


2  2 

>  8  sig  EV[s[x3]  s[x4]]  -  6  EV[s[x3]  ]  EV[s[x3]  s[x4]]  + 

2 

>  2  EV[s[x3]  s[x4]]  -  2  (2  EV[s[xl]  s  [x3]  ]  EV[s[xl]  s[x4]]  + 

2 

>  EV[s[xl]  ]  EV[s[x3]  s[x4]])  + 

2 

>  4  (EV[s[x3]  ]  EV[s[xl]  s[x4]]  +  2  EV[s[xl]  s[x3]]  EV[s[x3]  s[x4]])  + 

2  2  2  2  2  2 

>  2  sig  EV[s[x4]  ]  +  EV[s[xl]  ]  EV[s[x4]  ]  +  EV[s[x3]  ]  EV[s[x4]  ]  - 

2 

>  2  (2  EV[s[xl]  s[x4]]  EV[s[x3]  s[x4]]  +  EV[s[xl]  s[x3]]  EV[s[x4]  ]) 

4  2  2  2 

case3a=  6  sig  +  2  sig  EV[s[x2]  ]  +  4  sig  EV[s[x2]  s[x3]]  + 

2  2  2  2  2 

>  2  EV[s[x2]  s[x3]]  +  2  sig  EV[s[x3]  ]  +  EV[s[x2]  ]  EV[s[x3]  ]  - 

2  2 

>  8  sig  EV[s[x2]  s[x4]]  +  2  EV[s[x2]  s  [x4]  ] 

2  2 

>  8  sig  EV[s[x3]  s[x4]]  +  2  EV[s[x3]  s[x4]] 

2 

>  2  (2  EV[s[x2]  s[x3]]  EV[s[x2]  s  [x4]  ]  +  EV[s[x2]  ]  EV[s[x3]  s[x4]])  - 

2 

>  2  (EV[s[x3]  ]  EV[s[x2]  s[x4]]  +  2  EV[s[x2]  s  [x3]  ]  EV[s[x3]  s[x4]])  + 

2  2  2  2  2  2 

>  8  sig  EV[s[x4]  ]  +  EV[s[x2]  ]  EV[s[x4]  ]  +  EV[s[x3]  ]  EV[s[x4]  ]  - 

2  2 

>  6  EV[s[x2]  s[x4]]  EV[s[x4]  ]  -  6  EV[s[x3]  s[x4]]  EV[s[x4]  ]  + 

2  2 

>  3  EV[s[x4]  ]  +  4  (2  EV[s[x2]  s  [x4]  ]  EV[s[x3]  s  [x4]  ]  + 


2 


>  EV[s[x2]  s[x3]]  EV[s[x4]  ]) 

4  2  2  2 

case4a=  4  sig  +  2  sig  EV[s[xl]  ]  -  4  sig  EV[s[xl]  s[x2]]  + 

2  2  2  2 

>  2  sig  EV[s[x2]  ]  +  2  EV[s[xl]  s  [x3]  ]  +  2  EV[s[x2]  s  [x3]  ]  + 

2  2  2  2  2  2 

>  2  sig  EV[s[x3]  ]  +  EV[s[xl]  ]  EV[s[x3]  ]  +  EV[s[x2]  ]  EV[s[x3]  ]  - 

2 

>  2  (2  EV[s[xl]  s[x3]]  EV[s[x2]  s[x3]]  +  EV[s[xl]  s[x2]]  EV[s[x3]  ])  + 

2  2  2 

>  2  EV[s[xl]  s[x4]]  +  2  EV[s[x2]  s  [x4]  ]  -  4  sig  EV[s[x3]  s  [x4]  ]  - 

2 

>  2  (2  EV[s[xl]  s[x3]]  EV[s[xl]  s[x4]]  +  EV[s[xl]  ]  EV[s[x3]  s[x4]])  + 

>  4  (EV[s[x2]  s[x3]]  EV[s[xl]  s[x4]]  + 

>  EV[s[xl]  s[x3]]  EV[s[x2]  s[x4]]  +  EV[s[xl]  s[x2]]  EV[s[x3]  s[x4]])\ 

>  -  2  (2  EV[s[x2]  s[x3]]  EV[s[x2]  s  [x4]  ]  + 

2  2  2 

>  EV[s[x2]  ]  EV[s[x3]  s[x4]])  +  2  sig  EV[s[x4]  ]  + 

2  2  2  2 

>  EV[s[xl]  ]  EV[s[x4]  ]  +  EV[s[x2]  ]  EV[s[x4]  ]  - 

2 

>  2  (2  EV[s[xl]  s[x4]]  EV[s[x2]  s  [x4]  ]  +  EV[s[xl]  s  [x2]  ]  EV[s[x4]  ]) 

4  2  2  2 

caselb=  12  sig  +  24  sig  G[0]  +8  G[0]  -  24  sig  G[-x3  +  x4]  - 

2  2  2 

>  24  G[0]  G[-x3  +  x4]  +  4  G  [-x3  +  x4]  +4  (G[0]  +  2  G  [-x3  +  x4]  ) 

4  2  2  2 

case2b=  6  sig  +  12  sig  G[0]  +  6  G[0]  -  8  sig  G[-xl  +  x3]  - 

2  2 

>  6  G[0]  G[-xl  +  x3]  +  2  G[-xl  +  x3]  +  4  sig  G[-xl  +  x4]  + 


2 


2 


>  2  G[-xl  +  x4]  -  8  sig  G[-x3  +  x4]  -  6  G[0]  G[-x3  +  x4]  + 

2 

>  2  G[-x3  +  x4]  -  2  (2  G[-xl  +  x3]  G[-xl  +  x4]  +  G  [0]  G  [-x3  +  x4] )  + 

>  4  (G[0]  G[-xl  +  x4]  +  2  G[-xl  +  x3]  G  [-x3  +  x4] )  - 

>  2  (G[0]  G[-xl  +  x3]  +  2  G[-xl  +  x4]  G  [-x3  +  x4] ) 

4  2  2  2 

case3b=  6  sig  +  12  sig  G[0]  +  6  G[0]  +4  sig  G[-x2  +  x3]  + 

2  2 

>  2  G[-x2  +  x3]  -  8  sig  G[-x2  +  x4]  -  6  G  [0]  G[-x2  +  x4]  + 

2  2 

>  2  G[-x2  +  x4]  -  8  sig  G  [-x3  +  x4]  -  6  G [0]  G[-x3  +  x4]  + 

2 

>  2  G[-x3  +  x4]  -  2  (2  G  [-x2  +  x3]  G  [-x2  +  x4]  +  G [0]  G  [-x3  +  x4] )  - 

>  2  (G[0]  G[-x2  +  x4]  +  2  G  [-x2  +  x3]  G  [-x3  +  x4]  )  + 

>  4  (G[0]  G[-x2  +  x3]  +  2  G  [-x2  +  x4]  G  [-x3  +  x4]  ) 

4  2  2  2 

case4b=  4  sig  +  8  sig  G[0]  +  4  G[0]  -  4  sig  G[-xl  +  x2]  + 

2  2 

>  2  G[-xl  +  x3]  +2  G[-x2  +  x3]  - 

2 

>  2  (G[0]  G[-xl  +  x2]  +  2  G[-xl  +  x3]  G  [-x2  +  x3] )  +2  G[-xl  +  x4]  + 

2 

>  2  G[-x2  +  x4]  -  2  (G[0]  G[-xl  +  x2]  +  2  G[-xl  +  x4]  G  [-x2  +  x4] )  - 

2 

>  4  sig  G[-x3  +  x4]  -  2  (2  G[-xl  +  x3]  G[-xl  +  x4]  + 

>  G[0]  G[-x3  +  x4])  -  2  (2  G  [-x2  +  x3]  G  [-x2  +  x4]  + 

>  G[0]  G[-x3  +  x4])  +4  (G[-x2  +  x3]  G[-xl  +  x4]  + 

>  G[-xl  +  x3]  G [-x2  +  x4]  +  G[-xl  +  x2]  G [-x3  +  x4]  ) 

4  2  2  2 

caselc=  12  sig  +  24  sig  G[0]  +  12  G[0]  -  24  sig  G[p]  -  24  G[0]  G[p]  + 


2 


>  12  G[p] 

4  2  2  2 

case2c=  6  sig  +  12  sig  G[0]  +6  G[0]  -  16  sig  G[p] 

2  2 

>  12  G[p]  +4  sig  G[2  p]  +  4  G[0]  G [2  p]  -  8  G[p] 

4  2  2  2 

case3c=  6  sig  +  12  sig  G[0]  +6  G[0]  -  16  sig  G[p] 

2  2 

>  12  G[p]  +4  sig  G[2  p]  +  4  G[0]  G [2  p]  -  8  G[p] 

4  2  2  2 

case4c=  4  sig  +  8  sig  G[0]  +  4  G[0]  -  8  sig  G[p]  - 

2  2 

>  4  G[p]  +8  G[-xl  +  x3]  -  8  G[-xl  +  x3]  G[-x2  + 

2 

>  2  G  [-x2  +  x3]  -  8  G[-xl  +  x3]  G[-xl  +  x4]  + 

2 

>  4  G[-x2  +  x3]  G[-xl  +  x4]  +  2  G[-xl  +  x4] 


-  16  GEO]  G[p]  + 

2 

G[2  p]  +  2  G[2  p] 

-  16  GEO]  GEp]  + 

2 

GE2  p]  +  2  GE2  p] 

8  GEO]  GEp]  + 

x3]  + 


Appendix  D.  SSF  SNR  Results  for  DIMM  Geometry  H-  WFS 


N 


a  / 
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Atmospheric  turbulence  parameters,  such  as  Fried’s  coherence  diameter,  the  outer  scale  of  turbulence,  and  the 
turbulence  power  law,  are  related  to  the  wavefront  slope  structure  function  (SSF).  The  SSF  is  dehned  as  the 
second  moment  of  the  wavefront  slope  difference  as  a  function  of  both  time  and  position.  Knowledge  of  the 
SSF  allows  turbulence  parameters  to  be  estimated.  Hartmann  wavefront  sensor  (H-WFS)  slope  measurements, 
composed  of  both  signal  and  noise,  allow  the  SSF  to  be  estimated  by  computing  a  mean  square  difference  of 
H-WFS  slope  measurements.  The  quality  of  the  SSF  estimate  is  quantihed  by  the  signal-to-noise  ratio  (SNR)  of 
the  estimator.  This  thesis  develops  a  theoretical  SNR  expression  for  the  SSF  estimator.  This  SNR  is  a  function 
of  H-WFS  geometry,  the  number  of  temporal  frames  included  in  the  estimate,  the  outer  scale,  power  law,  and 
temporal  properties  of  the  turbulence.  Spatial  slope  correlations  are  incorporated.  Temporal  slope  correlations 
are  incorporated  using  Taylor’s  frozen  flow  hypothesis.  Results  are  presented  for  various  H-WFS  conhgurations 
and  atmospheric  turbulence  properties. 
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